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ABSTRACT 

We investigate the properties of the intracluster medium (ICM) that forms within 
-/V-body/hydrodynamical simulations of galaxy clusters in a ACDM cosmology. When 
radiative cooling and a simple model for galactic feedback are included, our clusters 
have X-ray luminosities and temperatures in good agreement with observed systems, 
demonstrating the required excess entropy in their cores. More generally, cooling and 
feedback increases the entropy of the ICM everywhere, albeit without significantly 
affecting the slope of the profile {S oc r) at large radii. The temperature of the ICM is 
only modestly increased by these processes, with projected temperature profiles being 
in reasonable agreement with the observations. Star/galaxy formation is still too effi- 
cient in our simulations, however, and so our gas mass fractions are around 60 per cent 
of the observed value at r25oo- Finally, we examine the reliability of using the hydro- 
static equilibrium equation to estimate cluster masses and find that it underpredicts 
the true mass of our clusters by up to 20 per cent, due to incomplete thermalisation 
of the gas. Feedback reduces this discrepancy, however, with estimates being accurate 
to within 10 per cent out to rsoo- 
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1 INTRODUCTION 

Clusters of galaxies play an important role in our quest for 
determining the cosmology of our Universe. In models of hi- 
erarchical structure formation such as the ACDM model, the 
mass function of clusters is a sensitive function of cosmolog- 
ical parameters, notably the matter density parameter, fim 
and the amplitude of linear fluctuations on 8 /i~^Mpc scales, 
erg. As a result, many authors have combined the observed 
(in the X-ray) abundance of clusters with theoretical mass 
functions (Press & Schechter 1974; Sheth, Mo & Tormen 
2001; Jenkins et al. 2001) to constrain these parameters (e.g. 
Evrard 1989; Henry & Arnaud 1991; Oukbir & Blanchard 
1992; White, Efstathiou & Frenk 1993; Viana & Liddle 1996; 
Eke, Cole & Frenk 1996; Henry 2000; Pierpaoli et al. 2003; 
Viana et al. 2003). 

In the post- WMAP era (Spergel et al. 2003; Verde 
2003), particular attention is now being paid to systematic 
uncertainties that may bias these estimates. As was recently 
pointed out by Henry (2004), the largest systematic uncer- 
tainty in determining as, using X-ray clusters is in the cali- 
bration of the X-ray temperature to mass relation. In order 
for robust determinations of such quantities to be made. 
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accurate theoretical predictions for how clusters form and 
evolve are required, taking full account of the physics that 
governs the structure of the intracluster medium (ICM). 

The best tool we have at our disposal to model the ICM 
is direct numerical simulation. Evrard (1990) performed the 
first cosmological simulation of a cluster and in his pio- 
neering study, demonstrated that many of its properties 
resembled those derived from X-ray observations. Subse- 
quent studies that followed (e.g. Thomas & Couchman 1992; 
Katz & White 1992; Cen & Ostriker 1994; Navarro, Frenk 
& White 1995; Evrard, Metzler & Navarro 1996; Bryan & 
Norman 1998; Eke, Navarro & Frenk 1998; Frenk et al. 1999; 
Thomas et al. 2001) have significantly improved our under- 
standing of X-ray clusters within hierarchical models. 

Many of the first cluster simulations ignored radia- 
tive cooling of the gas on the grounds that clusters today 
have average cooling times longer than the age of the uni- 
verse. These 'non-radiative' simulations vindicated simple 
self-similar scaling relations, as expected when the structure 
of the ICM is determined solely by gravity (Kaiser 1986). As 
useful as non-radiative simulations have been however, they 
contradict X-ray observations of clusters. The clearest exam- 
ple is the luminosity-temperature (Lx — Tx) relation, that 
is observed to have a steeper slope (Lx c>c Tx"^, e.g. Edge 
& Stewart 1991) than predicted from gravitational heat- 
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ing models (Lx oc Tx)- Another useful viewpoint is that 
the ICM is less concentrated because its entropy is higher 
(Kaiser 1991; Evrard & Henry 1991), as demonstrated by 
Ponman, Cannon & Navarro (1999). Consequently, much 
attention has been paid to models which focus on the en- 
tropy of the ICM (see e.g. Bower 1997; Voit & Bryan 2001; 
Voit et al. 2002,2003a). 

Early models concentrated on the preheating hypothe- 
sis, that the ICM was heated externally when the density 
contrast was low, therefore requiring less energy to reach 
the desired entropy level than if the heating occurred in a 
denser environment. Simulations of preheating models have 
been shown to be capable of producing the desired effect on 
the overall entropy level of the ICM (e.g. Navarro et al. 1995; 
Bialek, Evrard & Mohr 2001; Borgani et al. 2002) but have 
since fallen out of favour for many reasons, for example they 
predict isentropic cores in low-mass clusters that are not ob- 
served (sec Ponman, Sanderson & Finogucnov 2003). 

More recently, attention has shifted to modelling the 
eflfects of radiative cooling on the ICM. Naively, since the 
entropy of cooling gas decreases, one would expect cooling 
to have the opposite of the desired effect on the ICM, making 
clusters even more luminous than when cooling is neglected. 
Indeed, this effect was seen in some numerical studies (Katz 
& White 1992; Suginohara & Ostriker 1998; Lewis et al. 
2000), but was caused by the presence of an unrealistically- 
large central galaxy. Pearce et al. (2000) demonstrated that 
this 'overcooliug' problem could be ameliorated by 'decou- 
pling' the hot gas from the cold galactic gas, and as a result, 
found that cooling actually increased the ICM entropy. Most 
of the galactic material cools very quickly at high redshift 
without ever reaching the virial temperature (Binncy 1977; 
Kay et al. 2000; Binncy 2004), and instead emits in the 
Ultra- Violet (Fardal et al. 2000). The ICM adjusts itself to 
the loss of pressure support, with higher entropy gas flow- 
ing inwards to replace the cooled gas (Pearce et al. 2000). 
Various subsequent studies of cooling have successfully re- 
produced various observed X-ray properties of clusters (e.g. 
Bryan 2000; Muanwong et al. 2001,2002; Wu & Xue 2002; 
Dave, Katz & Weinberg 2002; Thomas et al. 2002; Voit et 
al. 2002). 

Although cooling contributes to the excess entropy of 
the ICM, heating processes (from stars and/or active galac- 
tic nuclei) are still required in order to regulate galaxy for- 
mation and possibly shape their luminosity function (e.g. 
White & Frenk 1991; Cole 1991; Balogh et al. 2001; Benson 
et al. 2003; Binney 2004). Numerical efforts are now focusing 
on incorporating such feedback processes with cooling into 
simulations of groups and clusters (Kay, Thomas & Theuns 
2003; Tornatore et al. 2003; Kay 2004; Borgani et al. 2004) 
and will continue to remain a vital area of study for the 
foreseeable future, until the effects of galaxy formation can 
be successfully incorporated into ICM models. 

In this paper we focus on a sot of high-resolution cos- 
mological simulations of 15 galaxy clusters, with the pri- 
mary aim of investigating how the combined effects of ra- 
diative cooling and a simple model for targeted heating by 
galactic feedback affects the structure of the ICM and as 
a consequence, estimates of the total mass distribution in 
clusters. Our study is in many ways similar to, and can be 
viewed as a progression of, those performed by Ascasibar et 
al. (2003) and Rasia, Tormen & Moscardini (2004), who re- 



stricted their simulations to non-radiative systems. We com- 
pare our findings with their results and with observational 
data where appropriate. 

The rest of this paper is organised as follows. In Sec- 
tion 2 we outline our method and discuss general properties 
of our simulated clusters. Radial profiles derived from the 
distribution of cluster mass (total mass density, ICM den- 
sity, entropy and gas fraction) are discussed in Section 3. 
In Section 4 we present cluster temperature profiles. Cluster 
mass determinations arc investigated under various condi- 
tions in Section 5. Finally, we draw conclusions in Section 6. 



2 SIMULATION DETAILS 

Simulations were performed assuming a ACDM cosmology, 
setting fin, = 0.3, Qa = 0.7, fib = 0.045, /i = 0.7 & era = 0.9. 

Our choice of parameters are in reasonably good agreement 
with the WMAP results (Spergel et al. 2003). 

All simulations were performed with version 2 of the 
GADGET A/^-body/hydrodynamics code (Springel, Yoshida & 
White 2001), kindly made available to us by V. Springel. 
Gravitational forces were computed using the combination 
of particle-mcsh and tree algorithms. Hydrodynamical forces 
were calculated using the variant of the Smoothed Particle 
Hydrodynamics (SPH) algorithm proposed by Springel & 
Hernquist (2002) that explicitly conserves entropy where the 
(artificial) viscosity of the flow is zero. 

Two sets of simulations were used in this study. First 
of all, we performed simulations of a random volume, sam- 
pling a cube of comoving length, L = 120/i~^Mpc us- 
ing N = 2 X 256^ particles (half baryons, half dark mat- 
ter). The dark matter and baryon particle masses are then 
mdark = 7.3 X lO'' h~^MQ and mbary = 1.3 x lO" /i'^Mq re- 
spectively. The second set of simulations consisted of 10 res- 
imulatod clusters, selected from a large A'"-body simulation 
performed by the Virgo Consortium (Yoshida, Slieth & Di- 
aferio 2001), with N = 512^ and L = 479 ft^^Mpc. The ten 
most massive objects with virial masses below 10^® h~^MQ 
were selected and resimulated at higher resolution, with sim- 
ilar particle masses to the random volume simulations. This 
sample has already been studied elsewhere without gas (see 
e.g. Gao et al. 2004). 

For all simulations, we started the runs at z = 49 and 
evolved them to 2 = (our results are presented for this 
redshift). We fixed the gravitational softening in comoving 
co-ordinates at all times to an equivalent Plummer value of 
e = 20 h^^kpc (gadget uses a spline softening kernel where 
the force becomes exactly Newtonian for r > 2.8e). 

2.1 The models 

In this paper we focus on two models, differing only in the 
freedom allowed for the entropy ^ of the gas to vary. In 
the first model, hereafter referred to as the Non-Radiative 
model, entropy can only increase due to shock-heating. Al- 
though it has already been demonstrated that this model 

^ We define entropy as S = kT (p/fimny ^, where 7 = 5/3 is 
the ratio of specific heats for a monatomic ideal gas and nmji = 
0.6 is the mean atomic weight of a fully-ionised plasma. 
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does not produce clusters that resemble observed systems 
(e.g. Evrard & Henry 1991) it nevertheless has been a well 
studied model and provides a good reference to measure the 
effects of non-gravitational processes on the ICM. 

In the second model, hereafter the Feedback model, we 
included two additional physical processes, both thought to 
be crucial in producing realistic clusters. First of all, the 
gas was allowed to lose entropy through radiative cooling. 
Although clusters at low redshift have average cooling times 
that arc longer than the age of the Universe, the same cannot 
be said for their precursors at high redshift when most of 
their galaxies form. The result of this removal of low entropy 
gas is a net increase in the entropy of tlu; remaining ICM 
(Pearcc et al. 2000). We included radiative cooling using 
the same method as used in the hydra code (see Thomas 
& Couchman 1992), adopting a tabulated cooling function 
(Sutherland & Dopita 1993) with Z = 0.32©. 

Gas is able to cool down to lO^K, where we arbitrarily 
cut off the cooling rate. The physics that governs the fate 
of this cooled gas, now identified as galactic material (e.g. 
Kay et al. 2000), is complex and would not be resolved in 
our simulations (although see e.g. Springel & Hernquist 2003 
for an attempt to model a multiphase galactic component 
within cosmological simulations). We instead adopt a phe- 
nomenological approach to galaxy formation, and assume 
that a fixed mass fraction of cooled gas will go on to form 
stars and the rest will be reheated by the stars (the heat- 
ing could also be due to active galactic nuclei, but note our 
model does not allow us to discriminate between the two). 

First, we identify cooled gas using a density thresh- 
old, nn > n* (where nn = Xp/mn is the hydrogen den- 
sity and X = 0.76 the hydrogen mass fraction) and a tem- 
perature threshold, T < T, , setting n, = 10^''' cm^'' and 
T* = 1.2x 10*K. Although typical densities where star/black 
hole formation occurs are much higher, we deliberately set 
the threshold to be low as it ensures that a negligible amount 
of thermal energy in the reheated gas is lost as radiation. 
Note therefore that we will be ignoring some cold gas that 
would otherwise be ablated by the ICM. 

Whether a particle then forms stars or is reheated is 
determined on a stochastic basis. The mass fraction of re- 
heated gas is controlled by a parameter /heat (hence a frac- 
tion, 1 — /heat will form stars) . For each cooled gas particle 
we draw a random number, r, from the unit interval and 
increase its entropy by a fixed amount, Shcat, if r < /hoat. 
The temperature of the reheated gas is then 

^^--=(100^) ( 10^)^^^ « 

We adopted a choice of parameters equivalent to that used 
by Kay (2004) in his Strong Feedback model, namely Sheat = 
lOOOkeVcm^ and /heat =0.1. For gas with density, rin = n*, 
'cT'heat = 17 keV, although part of this energy is distributed 
as the gas does work on its neighbours. 

This energy is only carried by 10 per cent of the mass 
of cooled gas, so it is useful to estimate the overall energy 
deposited into the ICM. Such an estimate is 

EiCM = Ti Afheat, (2) 

2 /urriH 

where Mheat = (/h^tt ~ l)~^Af« is the mass of reheated 
gas, directly related to the residual mass of stars in 



the cluster, M, . For our choice of parameters, -Bicm ~ 
1O®^(M./1O"M0) erg, comparable with the maximum en- 
ergy available from Type II supernovas, assuming one su- 
pernova releases 10®^ erg and the number of Type II super- 
novae, A^sNii = O.O1(M,/Af0). Thus, on average, the energy 
added to each ICM atomic nucleus is feTicM = 3.1(/*//gas) 
keV, where /» and /gas are the star and ICM mass frac- 
tions respectively. For clusters in our Feedback model, we 
find fcTicM ~ IkeV, i.e. a similar level to that found in pre- 
vious studies including heating (e.g. Wu, Fabian & Nulscn 
2000; Bower et al. 2001; Borgani et al. 2001; Muanwong et 
al. 2002). 



2.2 Cluster identification 

Clusters were identified using the procedure detailed in 
Muanwong et al. (2002). This involved finding the cluster 
centre (by searching for the density maximum) then growing 
spheres around this centre until the average density reached 
a constant factor. A, times the critical density. This thresh- 
old defines the relation between cluster mass and radius 

MA = ^Apcri?A, (3) 

where per = SHq/StyG is the critical density and Ho = 
70kms"^Mpc"^ for our cosmology. We define the virial 
mass/radius using A = 178 fi^*^ - 104 (Eke et al. 1998). 
Our main results will refer to virial quantities but occasion- 
ally we will use other values of A where appropriate. 

We only considered clusters with Mvir > 4.3 x 
lO'^* /i^'^Mq, corresponding to an effective particle number 
threshold, N = Mvir/(rn.dark + ?ngas) > 50,000, as suggested 
by Borgani et al. (2002) for obtaining reliable X-ray prop- 
erties. All our resimulated clusters and the 5 most massive 
objects in our random volume simulations satisfied this cri- 
terion, bringing our sample to 15 clusters per model. 



2.3 Cluster sample and global properties 

Table 1 lists global properties of our 15 Feedback clusters. 
The systems span a factor of 4 in mass, with fcTyir ~ 4 — 
SkeV. We define the ICM as all gas particles within R^ir with 
kT > O.lkeV (i.e. ~ 1.2 x lO^K). Generally, very little or no 
gas exists in our clusters below this temperature threshold 
(in the Feedback model, a few gas particles exist that have 
cooled down to T = 10'*K but have not yet formed stars or 
been reheated). 

There is very little scatter between ICM gas fraction 
values, /gas = 0.100 ± 0.005 (for the Non — Radiative clus- 
ters, /gas = 0.136±0.004). We note the average star fraction, 
/, = 0.035 ± 0.002, i.e. ~ 25 per cent of the baryon mass 
is coUisionless. By construction, this fraction is lower than 
would be expected in the absence of feedback, however it 
is still higher than the observed stellar mass fraction, mea- 
sured to be < 15 per cent (e.g. Balogh et al. 2001; Lin, Mohr 
& Stanford 2003). The observed value is an underestimate 
of /« as part of this mass may be hidden in a difi^use compo- 
nent (e.g. Murante et al. 2004), however note that the de- 
sired amount of cooled baryons can be achieved in our model 
by increasing /heat , without significantly affecting the global 
properties of the ICM (Kay et al. 2003). 
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Table 1. Global parameters of the simulated cluster sample, from 
left to right: cluster number; virial radius, in h~^Mpc; virial mass 
in 10^* /i~^Mq ; virial temperature in keV; number of baryon par- 
ticles (gas or stars) within iiviri ICM gas mass fraction; stellar 
mass fraction 



No. 


Rvir 




kT 




/gas 


/* 


1 


2.41 


16.83 


7.70 


184,082 


0.105 


0.035 


2 


2.03 


10.08 


4.91 


92,833 


0.102 


0.032 


3 


2.02 


9.82 


6.06 


88,840 


0.097 


0.035 


4 


1.99 


9.44 


5.57 


86,347 


0.098 


0.035 


5 


1.99 


9.43 


6.00 


89,872 


0.105 


0.034 


6 


1.96 


9.03 


5.10 


84,546 


0.103 


0.033 


7 


1.94 


8.68 


6.57 


82,922 


0.104 


0.035 


8 


1.82 


7.24 


4.34 


75,390 


0.098 


0.036 


9 


1.80 


6.93 


4.73 


68,042 


0.108 


0.035 


10 


1.79 


6.91 


4.66 


62,766 


0.097 


0.035 


11 


1.78 


6.80 


4.97 


64,532 


0.104 


0.034 


12 


1.73 


6.22 


4.02 


56,027 


0.096 


0.035 


13 


1.70 


5.88 


3.74 


62,316 


0.100 


0.036 


14 


1.67 


5.56 


4.54 


56,040 


0.091 


0.038 


15 


1.59 


4.77 


3.99 


48,156 


0.091 


0.038 




kT,/keV 

Figure 1. X-ray luminosity-temperature relation for clusters in 
the Non-Radiative (open squares) and Feedback (filled squares) 
models. Dashed lines are power law fits to the data (assuming 
Lx o<: for the Non-Radiative clusters). Symbols with error 
bars are for clusters studied by Markevitch (1998b) and Arnaud 
& Evrard (1998). 



To illustrate the difference between the X-ray proper- 
ties of the clusters in the two models, we show in Fig. 1 
the Lx — Tx relation at z = 0. Non-Radiative clusters are 
represented by open squares and Feedback clusters filled 
squares. X-ray emission- weighted temperatures and lumi- 
nosities were computed using the procedure outlined by 
Muanwong et al. (2002). To be approximately consistent 
with the observational data, luminosities are bolometric and 
temperatures calculated in a hard (1.5-8.0 keV) X-ray band. 
Furthermore, emission from within 50 /i~^kpc of the cluster 



□ □ A 



□ □ 

A A ° 
AAA 



kT,/keV 

Figure 2. Mass deposition rates for the Non- Radiative (squares) 
and Feedback (triangles) clusters, inferred from their core X-ray 
emission, as a function of emission- weighted temperature. Dashed 
lines give the mean rates. 



centre was omitted (i.e. we removed excess emission present 
due to cooling flows). 

The upper dashed line illustrates a best-fit power 
law to the Non-Radiative clusters, assuming 
This is the expected scaling relation for self-similar clus- 
ters assuming the dominant emission mechanism is thermal 
bremsstrahlung (Kaiser 1986); our clusters are in reasonable 
agreement with this scaling. 

Cooling and feedback processes combine to reduce 
the emission from within each cluster, particularly in low- 
temperature systems (see Kay 2004). This lowers and steep- 
ens the Z/x — 7x relation, bringing the Feedback clusters 
into good agreement with the observational data (Marke- 
vitch 1998b; Arnaud & Evrard 1998). The best-fit rela- 
tion for these systems is log(I/x,44) ~ ( — 1.42 ± 0.17) -I- 
(2.43 ± 0.23)log(ifcrx/keV), where Lx 44 is Lx in units of 
10"'' h-2ergs-\ 

Recently, it has emerged from X-ray spectra that there 
is a lack of cool gas (below one third of the mean tempera- 
ture) in the cores of clusters, even though the temperature of 
the gas is decreasing radially inwards and its cooling time is 
significantly shorter than a Hubble time (e.g. Fabian 2003). 
Such spectroscopic mass-deposition rates, Mspec, are of or- 
der 10 M0yr~'^, about a factor of 5-10 lower than rates 
estimated directly from the core X-ray emission of a cluster 

2 Lx fJ-mu 
5 



M: 



kTx 



(4) 



Unfortunately our model clusters do not have the resolution 
for us to study their cooling flows in detail (a lOMQyr"^ 
cooling flow corresponds to one gas particle cooling approx- 
imately every 200 Myr), in particular, to measure A/spcc- 



However, we can estimate the apparent mass deposition rate, 
Mx- Fig. 2 illustrates such rates for our Non-Radiative and 
Feedback clusters (calculated within 50 h~^'kpc of the cluster 
centre). Note that both models exhibit a large spread in val- 
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ues: in the former case, Mx ranges from 20 — 400 M0 yr~^ 
(although the gas can never cool), while in in the latter 
case, Mx = 3 — 13OM0yr~^. The Feedback model, with 
(Mx) =4OM0 yr ^, is in broad agreement with the obser- 
vations. 



3 CLUSTER MASS PROFILES 

Measuring the mass (both baryonic and total) distribution 
in clusters is of primary interest as it allows constraints to 
be placed on cosmological parameters, for example by cal- 
ibrating the relation between mass and X-ray temperature 
(e.g. Ifenry 2004) or luminosity (e.g. Allen et al. 2003), or 
by using the gas mass fraction itself (e.g. White et al. 1993b; 
Evrard 1997; Allen et al. 2002,2004). Knowledge of the ICM 
temperature distribution allows gas mass density profiles to 
be estimated via deconvolution of X-ray surface brightness 
profiles, and total masses estimated via the equation of hy- 
drostatic equilibrium. We will return to the latter issue in 
Section 5, but for now, examine the mass distribution of our 
simulated clusters directly. 



3.1 Mass density profiles 

Shown in Fig. 3 are total mass density profiles, plotted as 
x^p{x)/pcT, where pcr is the critical density and x — r/rvir- 
We set the minimum radius to be x^^in = 0.04 as this radius 
is both larger than the spline gravitational softening length 
and contains at least 100 hot gas particles for all clusters in 
our sample (Borgani et al. 2002; Ascasibar et al. 2003). Solid 
curves represent the mean profile and the error on the mean 
while error bars illustrate standard deviations within each 
bin, highlighting cluster to cluster variations. The scatter 
in density about the mean at each radius is quite small for 
both models, ~ ±15 per cent when averaged over all bins. 

A useful one-parameter model was proposed by Navarro 
et al. (1995, 1997, hereafter NFW) 



Non — Radiative 



p{x) _ Sc 

Per CX{1 + CXY ' 

where c is the concentration parameter and 



3 In(H-c) -c/(l-hc)' 



(5) 



(6) 



where A = 104 in our case. Best-fit NFW profiles to the 
mean density profiles, weighted by the error on the mean, 
are plotted as dashed curves in the figure. The mean profile 
concentration is 5.7 for the Non-Radiative model and 5.1 for 
the Feedback model. Cooling and feedback slightly fiatten 
the inner profile {x < 0.3), causing a ~ 10 per cent increase 
in the best-fit value of c. 

In the lower panels of Fig. 3, we plot the fractional dif- 
ference, Ap/p, between individual profiles and their corre- 
sponding NFW model fits. Again, the solid curve is the result 
for the mean profile and the error bars the standard devia- 
tion within each bin. Thomas et al. (2001) pointed out that, 
for their non-radiative clusters simulated in a rCDM cos- 
mology, the NFW model generally does not provide a good 
fit to individual systems but does well when fitting a profile 
averaged over clusters with similar mass. For our clusters. 



1 

0.2 
0.0 
< -0.2 




c = 5.i 



Q. 



a = 0^028 ^ 
Feedback 



Q. 



Q. 




c = 5.12 



.2 r 

.0 

.2 r 



a = 0.040 



0. 



.0 



Figure 3. Total mass density profiles, x^p{x), where x = r/rvir, 
for clusters in the Non- Radiative and Feedback models. Solid 
curves are the mean profile and the error on the mean, and 
the error bars illustrate the standard deviation within each bin. 
The dashed curve is the best-fit NFW profile to the mean curve, 
weighted by the error on the mean, with the concentration pa- 
rameter given in the legend. Also plotted below each main panel 
is the fractional diff'erences between individual profiles and cor- 
responding best-fit NFW profiles. Again, the solid curve is the 
mean and the error bars the standard deviations within each bin. 
The rms fractional difference, averaged over all bins, is given in 
the legend. 



the fit also improves significantly when averaging over clus- 
ters. The mean deviation at a given radius is around 10 
per cent, but the root-mean-square difference, {(Ap/p)^), 
when averaged over all clusters and bins is only a few per 
cent. 



3.2 ICM density profiles 

Fig. 4 illustrates ICM density profiles: comparing the two 
simulation models reveals that the Feedback clusters have a 
lower density than the Non-Radiative clusters everywhere, 
with the difference being particularly pronounced in the cen- 
tral region. This is caused both by cooling (which removes 
gas from the ICM) and feedback (which heats the ICM). 
Since the cooling rate scales as p and the feedback rate effec- 
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Non — Radiative 



10 7 




x=Om, (3 = 0.84 



Q. 



0.2 
0.0 



o -0.2 




a. .=-0.052 ^ 

Feedback 



10 7 



Q. 




x. = 0.10, (3 = 0.75 



.2 7 



.2 7 



a = 0.072 



0.1 



.0 



Figure 4. ICM density profiles, x'^pgss{x), fitted using the 
isotliermal /3-modeI. The model profile fitted to the results of 
Rasia et al. (2004) is shown as a dotted curve. 



tively scales with the cooling rate, we expect their influence 
to be more pronounced closer to the cluster centre. 

Historically, the isothermal /3-model (Cavaliere & 
Fusco-Femiano 1976) is used to fit X-ray surface bright- 
ness profiles (we will return to this in Section 5), with the 
underlying ICM density profile being 

Pgas(a;) _ <5/3 



fhPc 



(a:2-ha:i)3/'/2' 



(7) 



where /b — Sib /Jim, = rc/r^ii is the core radius and /3 
determines the asymptotic slope at large radii. We confirm 
in Fig. 4, for both our simulation models, the finding by 
other authors (e.g. Ascasibar et al. 2003, hereafter AYMG; 
Rasia et al. 2004, hereafter RTM) that the model is not a 
good fit to the data, systematically over-predicting the slope 
of the profile at small radii and under-predicting it at large 
radii. 

RTM suggested an alternative expression for the density 
profile, which we write as 



/bPcr 



(8) 



where they fixed Xp — 0.04 and a = 2.5 (shown in Fig. 4 as a 
dotted curve). As can be seen however, this profile does not 
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Figure 5. Dimensionless entropy profiles for clusters in the Non- 
Radiative (top) and Feedback (bottom) models. The dashed line 
is the best-fitting power law to the mean profile for 0.2 < x < 1 
(the power-law index is given in the legend). 



provide a good fit to our data. A good fit can be achieved 
for our Non-Radiative clusters if we instead allow the data 
to select values for Xp and a; our results favour a steeper 
asymptotic slope [a ^ 3) and the ICM in our clusters is less 
centrally concentrated than found by RTM. For the Feedback 
clusters, equation 8 also improves the fit over the /3-model, 
but still does not provide a satisfactory description of the 
data, again getting the slope wrong at large and small radii. 
We therefore caution the use of such expressions as a general 
tool in the modelling of clusters: a larger sample of clusters 
spanning a wider range in mass will be required in order to 
investigate this issue further. 



3.3 ICM entropy profiles 

The ICM entropy profile is a particularly useful quantity as 
it reveals information about the nature of non-gravitational 
processes in clusters (e.g. Evrard & Henry 1991; Kaiser 1991; 
Bower 1997; Voit & Bryan 2001; Voit et al. 2002,2003a; Pon- 
man et al. 2003; Voit & Ponman 2003b; Kay 2004; Borgani et 
al. 2004). Gravitational heating models, in which the outer 
entropy profile is determined by a spherical accretion shock, 
predicts an outer slope of 1.1 (Tozzi & Norman 2001). Ob- 
served entropy profiles exhibit similar values (e.g. Arnaud, 
Pratt & Pointecouteau 2004), suggesting that the main ef- 
fect of non-gravitational processes is to increase the normal- 
isation of the entropy profile without significantly affecting 
its shape (see Ponman et al. 2003). 

In Fig. 5 we present ICM entropy profiles, scaled by 
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Figure 6. ICM gas fraction (left) and baryon fraction (right) profiles for clusters in the Non- Radiative and Feedback models. The global 
baryon fraction, r2t,/nm = 0.15, is plotted as a horizontal dashed line. The triangle is the error-weighted mean measurement at r200 
{x ~ 0.75) by Ettori (2003). Square symbols are data points from Allen et al. (2004) 



5'vir = fcTvir(Mf"H)^^^(Apcrnb/nm)"^/^ The profiles are 
very close to power laws for x > 0.2. Our Non-Radiative 
clusters predict an outer slope of 1.2, very similar to that 
predicted by Tozzi & Norman's models, and consistent with 
a study of smaller systems, shown to be in good agreement 
with results from independent simulations using enzo, an 
adaptive mesh refinement code (Voit, Kay & Bryan 2004). 
Cooling and feedback increase the entropy distribution of 
the gas at all radii, acting to flatten the profiles: at x = 1 
the increase is only 25 per cent but rises to almost a factor of 
2 at a; = 0.2. Note, however, the outer slope does not change 
signiflcantly, S (x x, and so is still similar to observational 
determinations. 



3.4 Gas fraction profiles 

A unique feature of cluster-sized haloes is that, in principle, 
most (if not all) of their baryonic content can be observed di- 
rectly. Combined with an estimate of the total cluster mass, 
it is then possible to place constraints (in particular, an up- 
per limit) on f2m, assuming a baryon density inferred from 
elsewhere, such as the nucleosynthetic value (White et al. 
1993b). 

Panels on the left of Fig. 6 illustrate cumulative ICM 
gas fraction profiles, /gas(< x) = Mgas(< x)/Mtot{< x). As 
was found by previous authors (e.g. Eke et al. 1998; Frenk et 
al. 1999) the Non-Radiative profile appears to converge to a 
value close to, but not exactly, the global value at the virial 
radius (we find /gas ~ 0.9r2b/f^m)- (The baryon fraction does 



not reach the global value until x ~ 2.) The rise is gradual: 
/ga4a; = 0.1)~0.8/gaB(a; = l). 

Gas fractions in the Feedback clusters are 25 per cent 
lower at X = 1 due to radiative cooling: as can be seen from 
the right half of the figure, the total baryon fraction profiles 
reach the same value as in the Non-Radiative clusters. The 
increase in gas fraction with radius is also larger, /gas(a; — 
0.1) ~ 0.6/gas(a; = 1), as a result of the gas being more 
extended due to its higher entropy. 

Recently, Allen et al. (2004) measured gas fraction pro- 
files for 26 relaxed X-ray luminous clusters observed with 
Chandra, out to A = 2500. They found that the profiles 
rise from /gas ~ 0.07 at a; = to /gas ~ 0.12 at x = 0.25 
('"2500) • These results are also shown in Fig. 6, plotted as 
squares. Intriguingly, the observational data resemble our 
Non-Radiative clusters more than our Feedback clusters, 
with the latter having gas fractions that are around 60 per 
cent of the observed value at r25oo- We also show the error- 
weighted mean gas fraction at r2oo {x ~ 0.75), as measured 
by Ettori (2003), consistent with the result of Allen et al. 
if the profiles do not continue to rise much beyond r25oo- 
Here, our Feedback results are approximately 80 per cent of 
the observed estimate. 

In the right-hand panels of Fig. 6 we plot bias factor pro- 
files, b = /bary/(^^b/f2m), i.e. the baryon fraction, /bary, in 
units of the global value. Allen et al. assume a constant frac- 
tion (16 per cent) of baryons in stars. Interestingly, the Feed- 
back clusters are now in good agreement with their observa- 
tions. Cooling leads to a rise in the central baryon fraction, 
fiattening the profile with respective to the Non-Radiative 



© 2002 RAS, MNRAS 000, 1-15 



8 S. T. Kay et al. 



case. Clearly however, our Feedback clusters have too much 
material in stars at all radii (/* decreases monotonically 
with radius, from 0.54 at 2: = 0.05 to 0.25 at a: = 1). At 
least part of the problem may be due to forming stars at 
artificially low densities: while increasing the star formation 
density threshold presents other problems (in particular, the 
propagation of feedback energy to lower density regions), it 
will allow more gas to be ablated as it moves through the 
ICM, therefore reducing the amount of material available to 
form stars. 



4 TEMPERATURE PROFILES 

A more accurate description of the ICM requires the spatial 
distribution of its temperature to be measured, which can 
be achieved using X-ray spectroscopy. Both observations and 
simulations concur that the ICM is non-isothermal, and so 
measuring temperature gradients improves the accuracy of 
estimating a cluster's mass if it is in hydrostatic equilibrium. 

Early observations of azimuthally-averaged tempera- 
ture profiles debated whether the inner profile was declining 
with radius (Markevitch et al. 1998a) or was isothermal (Ir- 
win, Bregman & Evrard 1999; White 2000; Irwin & Breg- 
man 2000). A later study by De Grandi & Molendi (2002) 
was consistent with isothermality outside cooling-fiow re- 
gions and within ~ 0.2rvir. This was confirmed by Allen et 
al. (2001) for their sample of 6 luminous relaxed clusters 
observed with Chandra, and more recently by Arnaud et al. 
(2004) for 7 clusters observed with XMM-Newton. Although 
the situation is still not robust, the general consensus is that 
the ICM temperature rises from the centre to ~ O.lrvir, 
is approximately isothermal then starts to decline beyond 
~ 0.2 - O.Srvir. More high-quality data from XMM-Newton 
and Chandra will improve the situation considerably. 



4.1 3D temperature profiles 

We first present 3D temperature profiles, shown in Fig. 7. 
Each profile is scaled by the isothermal model virial tem- 
perature of the cluster, Tvir ~ GfimiiMviT /2kRvir. For the 
Non-Radiative clusters, our results are in good agreement 
with previous numerical studies (e.g. Evrard 1990; Navarro 
et al. 1995; Eke et al. 1998), i.e. the profile slowly varies 
out to 2; ~ 0.2, then declines by around a factor of 2 out to 
a; = 1. A more contentious issue is whether the inner profile 
is isothermal (e.g. RTM) or continues to rise to the cen- 
tre (e.g. AYMG) . We find our average Non-Radiative profile 
is isothermal within x ~ 0.1, although the core tempera- 
ture is hotter than that found by RTM (their best-fit profile 
is shown as a dotted curve). Given AYMG also used the 
entropy-conserving version of GADGET (RTM did not), the 
difference in the slope of the inner profile could be due to the 
selection of objects (AYMG studied lower-mass systems), 
although we note both RTM and AYMG performed studies 
with both higher mass and force resolution than ours. 

Including cooling and feedback produces a modest in- 
crease (~ 10 per cent at a; = 1) in the temperature of the 
ICM. Note also that the average Feedback profile begins to 
turn over within a; 0.1 due to radiative cooling. 




Figure 7. Dimensionless temperature profiles for clusters in the 
Non-Radiative and Feedback models. The best-fitting profile to 
the results of Rasia et al. (2004) is shown as a dotted curve. 



4.2 Projected temperature profiles 

To compare our temperature profiles with observations, 
we have constructed projected profiles. Here, we emulate 
the method employed by Loken et al. (2002), who them- 
selves tried to mimic the procedure employed by Marke- 
vitch et al. (1998a). X-ray emission- weighted temperature 
and surface brightness maps were first constructed for the 
1.5 — llkeV band, using the procedure detailed by Onuora, 
Kay & Thomas (2003). Each map was constructed from a 
cube of width 4/i~^Mpc, centred on the cluster, and the 
length of the pixels were set to the gravitational softening, 
e = 0.02 /i"^Mpc (hence each map is 200 x 200 pixels) . Maps 
were then centred on the pixel containing the maximum 
surface brightness and normalised by the X-ray emission- 
weighted temperature (Tx) within 1 h~^Mpc of this centre. 

The temperature maps (with surface brightness con- 
tours overlaid) are shown in Fig. 8, out to one Abell radius 
(1.5 /I'^Mpc). It is clear from this figure (see also Loken 
et al. 2002; Onuora et al. 2003) that the temperature of the 
ICM is very asymmetric and relates to the individual merger 
history of each cluster. How much of this information can 
be recovered from observed clusters, however, is a question 
that is now only starting to be addressed (e.g. Gardini et al. 
2004; Mazzotta et al. 2004). 

We then constructed projected emission-weighted 
temperature profiles by azimuthally averaging temper- 
atures within bins of width 0.05-Ri8o, where iiiso ~ 
1.95(fcrx/10keV)^/^ /i"^Mpc is the radius commonly used 
by observers, calibrated from the simulations of Evrard et 
al. (1996). The results are shown in the top panels of Fig. 9, 
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Figure 8. X-ray temperature maps for the clusters studied in this paper. Maps are shown within an Abell Radius (1.5 h ^Mpc). X-ray 
surface brightness contours (with order of magnitude intervals) are overlaid. 



compared to the data of Markevitch et al. (1998a); De 
Grandi & Molendi (2002) and Arnaud et al. (2004). Intrigu- 
ingly, our mean profiles are in reasonably good agreement 
with the observations within 0.2i?i8o (excluding the cooling- 
flow region in the Non-Radiative model) although they fall 
off less rapidly in the outer parts. This is contrary to recent 
claims by several groups (e.g. Loken et al. 2002; AYMG; 
Borgani et al. 2004), who found consistency with the obser- 
vations at large radii but whose profiles continued to rise 
into the centre. For example, Borgani et al., who included 
a model for both cooling and feedback in their calculation, 
found that the average temperature of clusters above 3keV 
increased to 1.2 times the emission-weighted temperature 
at ~ 0.3i?i8o, before turning over. It is unclear whether the 
discrepancy is due to numerical resolution effects (Borgani 
et al. simulated their clusters at higher resolution) or differ- 
ences in the feedback model or both. It is imperative that 



such discrepancies are fully understood and we aim to do so 
in future work. 

In the bottom panels of Fig. 9, we also compare pro- 
jected mass-weighted profiles to the result found by Allen 
et al. (2001), plotted within i?2500- Our Feedback model is 
in reasonably good agreement at these radii (to within ~ 10 
per cent). 



5 MASS ESTIMATES 

We now investigate the accuracy with which our cluster 
masses can be recovered when making various assumptions 
regarding the radial structure of the ICM. If the ICM is in 
hydrostatic equilibrium, the mass can be calculated from 



M(< x) 



Gp{x) dr 



(9) 
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Figure 9. Projected ICM temperature profiles compared with various observational data. Left panels: profiles are plotted out to 0.7i?i80 
and compared with results from Markevitch et al. (1998a; dashed box region), De Grandi & Molendi (2002; open and filled squares) and 
Arnaud, Pratt & Pointecouteau (2004; solid rectangle). Right panels: profiles are plotted out to i?2500 a^nd compared to the result from 
Allen et al. (2001). 



The simplest model would then be to assume the ICM is 
a polytropic gas, P = up'' . For isothermal distributions, 
7 = 1, and for adiabatic distributions, 7 = 5/3. We show 
in Fig. 10 (plotting the effective polytropic index, 7 = 1 + 
dln(r)/dln(p), at each radius), that this assumption is not 
valid in general for our simulated clusters. For the Non- 
Radiative systems, 7 ~ 1.1 out to a; ~ 0.2 then increases 
to 7 ~ 1.3 at a: = 1. For the Feedback systems, 7 increases 
monotonically from the centre outwards, varying from ~ 0.8 
at a; = 0.04 to ~ 1.3 at a; = 1. Note in this case, 7 is less 
than unity in the cluster core, reflecting the positive radial 
temperature gradient induced by cooling. 

AYMG split their sample of poor clusters into relaxed, 
minor and major merger systems and found the first two 
to be well described by a single polytrope with 7 = 1.18 
(shown in our figure as a horizontal dashed line). For their 
major mergers, the main deviation in 7 was in the central 
regions [x < 0.3), where the gas was close to isothermal. On 
average, our Non-Radiative clusters also have lower values 
of 7 in their centres, suggesting that our objects are less 
relaxed than theirs, as would be expected given their higher 
mass and hence longer dynamical times. 



5.1 The isothermal /3-model 

Although our clusters are clearly not isothermal, it is nev- 
ertheless interesting to apply the isothermal /3-model (Sec- 
tion 3.2) to equation 9 as it is this model that is most com- 



monly used by observers to estimate cluster masses. First, 
the X-ray surface brightness profile is fitted 

where Rc is the core radius and /3flt is the asymptotic slope 
of the profile at large radii. The results from applying the 
same procedure to our surface brightness maps are shown 
in Fig. 11. As has already been shown when fitting the ICM 
density profile (Fig. 4), it is clear that this model does not 
describe the simulation data well over the whole range of 
radii. We only fit the data for x < 0.5 (at larger radii the 
profiles diverge due to infalling substructure) ; even then, the 
average fractional difference is around 10 per cent. 

The total mass within a given radius can then be esti- 
mated 

where Tx is a single (emission-weighted) measurement of 
the ICM temperature. In Fig. 12 we show results from ap- 
plying equation 11 to our clusters, plotting the ratio of the 
estimated to the true mass within a given radius. The re- 
sults are very similar for the Non-Radiative and Feedback 
clusters: equation 11 underestimates the mass by ~ 20 per 
cent at a; < 0.3 and overestimates the mass by ~ 30 per 
cent at a: = 1. Interestingly, at a: = 0.5 (~ -R500) the esti- 
mated mass is within 5 per cent of the true mass. Within 
^^200 (a^ 0.75), masses can be estimated to better than 
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Figure 10. Effective polytropic index, 7, as a function of ra- 
dius for the Non- Radiative and Feedback models. Tlie horizontal 
dashed line marks the value 7 = 1.18, found to be a good de- 
scription for the relaxed poor clusters studied by Ascasibar et al. 
(2003). 



20 per cent. A similar study was performed by Evrard et 
al. (1996), who considered various models and found a sim- 
ilar trend with radius, albeit better agreement between the 
estimated and true masses at x < 0.5. 



5.2 Single versus spatial temperature models 

In the new era of Chandra and XMM-Newton, obtaining 
spatially resolved density and temperature information is 
now becoming possible. It is therefore of interest to estimate 
cluster masses when accurate descriptions for their density 
and temperature profiles are known. 

Following RTM, we first perform the test where we use 
a single temperature to describe the ICM but use its true 
density profile. Equation 9 becomes 



kT{< x)xrvir d\np{x) 
Gfirrm dln{x) ' 



(12) 



where r(< x) is the mass-weighted temperature within x. 
Results are shown in Fig. 13, again plotting the estimated 
to true mass ratio as a function of radius. 

Our Non-Radiative result is in reasonably good agree- 
ment with the findings of RTM, with the estimated mass 
being lower, but within 20 per cent of the true mass for 
0.1 < X < 1. Ata;<0.3, ratios for the Feedback clusters are 
similar, but the estimated mass approaches the true mass 
much more rapidly at larger radius, being within ~ 5 per 
cent between x — 0.5 (i?5oo) and x = 1. The integrated 




0.1 



Figure 11. X-ray surface brightness profiles for clusters in the 
Non-Radiative and Feedback models. Fits are only produced for 

R < O.SHvir. 



mass-weighted temperature is ~ 10 per cent larger in the 
Feedback clusters at all radii, increasing the estimated mass 
by the same factor. This is cancelled out at small radii by 
the flatter inner slope of the density profile. At larger radii 
however, the slope increases with radius, eventually approx- 
imating the same slope as in the Non-Radiative clusters and 
causing the mass estimate to rise with radius. 

We now include spatial temperature information in our 
mass estimate 



kT{x) X rv 
G/irriH 



dlnp dluT 



d\nx d In X 



(13) 



with results shown in the middle panels of Fig. 13. Our Non- 
Radiative results are again in good agreement with RTM. 
Both models underpredict the mass at all radii by up to 
20 per cent. Again, the mass estimate is slightly higher for 
the Feedback clusters because of the increase in ICM tem- 
perature. Clearly from these results, our clusters cannot be 
completely in hydrostatic equilibrium. As pointed out by 
RTM, the isothermal mass estimator is more accurate than 
the non-isothermal case at larger radii, because the inte- 
grated temperature happens to cancel out the effects of gas 
motion. 

The lack of agreement between estimated and true clus- 
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ter masses is due to the presence of kinetic support, which 
becomes particularly significant at large radii. Radial veloc- 
ity profiles confirm that our clusters are, to very good ap- 
proximation, virialised systems. However, as Fig. 14 shows, 
the velocity dispersion of the ICM is non- negligible and be- 
comes increasingly significant with radius. Note, however 
that kinetic support is less important for the Feedback clus- 
ters: the Non-Radiative clusters have a positive velocity dis- 



persion gradient at all radii, whereas in the Feedback clusters, 
a significant gradient appears at a; > 0.2. Thus, the ICM in 
the Feedback clusters is hydrostatic to better approximation 
than in the Non-Radiative clusters. 

Including an isotropic velocity dispersion term in equa- 
tion 9, namely P — p{kT / fimn + cr^/3), we show ratios of 
estimated to true mass estimates in the bottom panels of 
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Fig. 13. On average, the estimated masses are now in good 
agreement with the true masses. 



6 CONCLUSIONS 

In this paper we have studied high-resolution A'^- 
body/hydrodynamical simulations of galaxy clusters to in- 
vestigate the effects of radiative cooling and galactic feed- 
back on the properties of the ICM. In particular, we com- 
pared a set of Non-Radiative simulations (the simplest model 
that is well-studied in the literature, but is nevertheless ruled 
out by X-ray observations of clusters) to our Feedback model 
where a fraction (10 per cent) of cooled material is given a 
fixed amount of entropy (lOOOkeVcm^) that is then redis- 
tributed in the ICM. The latter model reproduces the ob- 
served cluster Lx — 2x relation, signifying the correct level 
of excess entropy in the ICM core. Our main conclusions are 
as follows. 

• The distribution of total mass is, on average, reasonably 
well described by an NFW profile, with residual differences 
within 5 per cent. Including cooling and feedback does not 
significantly alter the shape of the cluster potential, causing 



only a ~ 10 per cent increase in the best-fit concentration 
parameter. 

• Cooling and feedback act to increase the entropy of the 
ICM everywhere, although the effect is more pronounced 
within the central region, where the cooling rate (and hence 
the rate of mass/energy injection) is higher. Outside ~ 
0.2r-vir, the entropy profile scales approximately linearly with 
radius, consistent with recent X-ray observations. 

• As a consequence of the excess entropy, the ICM den- 
sity profile is lower and fiatter in the Feedback clusters. The 
isothermal /3-model does not provide a good fit to the ICM 
density profile in either numerical models, systematically 
overpredicting the inner slope and underpredicting the outer 
slope. 

• Baryon fractions rise gradually from the centre out- 
wards, reaching ~ 90 per cent of the global value (f2b/f2m) 
at the virial radius. Virtually no gas is lost from the Feedback 
clusters due to galactic outfiows; the lower ICM gas fraction 
can be accounted for by the fraction of stars. However, the 
Feedback model does not contain enough gas within r25oo 
to agree with the observed determinations by Allen et al. 
(2004), although the overall baryon fraction profiles are in 
good agreement. Our model overproduces the mass in stars 
at all radii, possibly a symptom of having such a low density 
threshold for star formation. 
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• Cooling and feedback produce only a modest increase in 
the temperature of the ICM (around 10 per cent at the virial 
radius) and cooling induces a positive radial temperature 
gradient for r < O.lrvir- The temperature of the ICM is then 
primarily driven by gravity. Projected temperature profiles 
are in reasonable agreement with the observations although 
the decline at larger radii is not as severe. 

• On average, the ICM cannot be described by a single 
polytropic equation of state: the effective polytropic index, 
7 = 1 + d\n{T)/d\n{p) increases with radius, reaching 7 ~ 
1.3 at the virial radius. Central values are lower than unity in 
the Feedback clusters, due to the positive radial temperature 
gradient induced by radiative cooling. 

• Estimating the total mass of clusters by fitting surface 
brightness profiles with the isothermal /3-model is accurate 
to within 20-30 per cent (depending on the outer radius 
chosen), regardless of the simulated model. The model is 
most accurate around r = O.Srvir (~ rsoo), where the error 
is ~ 5 per cent. 

• As found by Rasia et al. (2004), the ICM is not com- 
pletely hydrostatic, leading to an underestimate in the true 
mass of the cluster of up to 20 per cent in our Non-Radiative 
clusters, even when full density and temperature information 
is included in the calculation. Including a kinetic pressure 
term recovers the true mass. The degree of thermalisation in 
the ICM improves when cooling and feedback are included, 
predicting a hydrostatic mass that is within 10 per cent of 
the true mass for r < r^oo- 

In summary, we find that cooling and feedback mainly 
aS'cct the inner structure of the ICM, as a consequence of the 
excess entropy produced by these processes. Although our 
Feedback model approximately achieves the overall desired 
entropy excess in clusters, it does not provide a detailed 
match to the inner structure of the ICM, namely the gas 
mass fraction within r25oo is around 60 per cent of the ob- 
served value (Allen et al. 2004). Whether this discrepancy 
can be rectified within the context of our feedback model 
(by varying the model parameters and/or increasing the nu- 
merical resolution of the simulations), or whether our simu- 
lations axe still missing some fundamental physical process, 
requires further investigation. 
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